clear all

use "Estimates_for_plots_label.dta"

destring month, replace
*rename estimate coef


graph set window fontface "Times New Roman"

rename estimate coef

**HUNGARY, 2016
*++++++++++++++++++++

***Subsidized unit sales by labels A+++, A++, and A+, time path, window around reform [-3,1]. Longer window impossible, given data.
twoway(rarea max95 min95 month if type=="s" & program=="hu16" & label=="a3", color(gray%15)) (rarea max95 min95 month if type=="s" & program=="hu16" & label=="a2", color(green%50) lcolor(green%50) lpattern(dash)) (rarea max95 min95 month if type=="s" & program=="hu16" & label=="a1", color(none) lcolor(black) lpattern(dot)) (scatter coef month if type=="s" & program=="hu16" & label=="a3", connect(l) tline(0, lstyle(refline) lcolor(black) lwidth(0.1)) tline(3, lstyle(refline) lwidth(0.1) lcolor(black) lpattern(dash)) yline(0, lstyle(refline) lcolor(gs8) lwidth(0.1)) lpattern(solid) lcolor(black) mcolor(black) msym(D) msize(small) graphregion(color(white)) legend(ring(0) position(10) row(3) lab(4 "A+++")lab(1 "95% CI")lab(5 "A++") lab(2 "95% CI") lab(6 "A+")lab(3 "95% CI") region(lstyle(none) color(none)) nobox) leg(on order(4 1 5 2 6 3)) xtitle("Month") ytitle("log(UNITS)") ylabel(-1(1.5)2.5)) (scatter coef month if type=="s" & program=="hu16" & label=="a2", connect(l) lpattern(solid) lcolor(red) msize(small) mcolor(red) msym(O) xlabel(-3 "Jun."-2 "Jul." -1 "Aug." 0"Sept." 1 "Oct." 2 "Nov." 3"Dec." 4 "Jan.")) (scatter coef month if type=="s" & program=="hu16" & label=="a1", connect(l) msymbol(T) msize(small))

graph export "HU_16_label2.eps", as(eps)  fontface("Times New Roman") replace


***Subsidized prices by labels A+++, A++, and A+, time path, window around reform [-3,1]. Longer window impossible, given data.
twoway(rarea max95 min95 month if type=="sp" & program=="hu16" & label=="a3", color(gray%15)) (rarea max95 min95 month if type=="sp" & program=="hu16" & label=="a2", color(green%50) lcolor(green%50) lpattern(dash)) (rarea max95 min95 month if type=="sp" & program=="hu16" & label=="a1", color(none) lcolor(black) lpattern(dot)) (scatter coef month if type=="sp" & program=="hu16" & label=="a3", connect(l) tline(0, lstyle(refline) lcolor(black) lwidth(0.1)) tline(3, lstyle(refline) lwidth(0.1) lcolor(black) lpattern(dash)) yline(0, lstyle(refline) lcolor(gs8) lwidth(0.1)) lpattern(solid) lcolor(black) mcolor(black) msym(D) msize(small) graphregion(color(white)) legend(ring(0) position(10) row(3) lab(4 "A+++")lab(1 "95% CI")lab(5 "A++") lab(2 "95% CI") lab(6 "A+")lab(3 "95% CI") region(lstyle(none) color(none)) nobox) leg(on order(4 1 5 2 6 3)) xtitle("Month") ytitle("log(PRICE)") ylabel(-0.05(.05)0.1)) (scatter coef month if type=="sp" & program=="hu16" & label=="a2", connect(l) lpattern(solid) lcolor(red) msize(small) mcolor(red) msym(O) xlabel(-3 "Jun."-2 "Jul." -1 "Aug." 0"Sept." 1 "Oct." 2 "Nov." 3"Dec." 4 "Jan.")) (scatter coef month if type=="sp" & program=="hu16" & label=="a1", connect(l) msymbol(T) msize(small))

graph export "HU_16_label_p2.eps", as(eps)  fontface("Times New Roman") replace



**HUNGARY, 2015
*++++++++++++++++++++


***Subsidized unit sales by labels A+++, A++, and A+, time path, window around reform [-3,3].
twoway(rarea max95 min95 month if type=="s" & program=="hu15" & label=="a3", color(gray%15)) (rarea max95 min95 month if type=="s" & program=="hu15" & label=="a2", color(green%50) lcolor(green%50) lpattern(dash)) (rarea max95 min95 month if type=="s" & program=="hu15" & label=="a1", color(none) lcolor(black) lpattern(dot)) (scatter coef month if type=="s" & program=="hu15" & label=="a3", connect(l) tline(0, lstyle(refline) lcolor(black) lwidth(0.1)) tline(3, lstyle(refline) lwidth(0.1) lcolor(black) lpattern(dash)) yline(0, lstyle(refline) lcolor(gs8) lwidth(0.1)) lpattern(solid) lcolor(black) mcolor(black) msym(D) msize(small) graphregion(color(white)) legend(ring(0) position(10) row(3)  lab(4 "A+++")lab(1 "95% CI")lab(5 "A++") lab(2 "95% CI") lab(6 "A+")lab(3 "95% CI") region(lstyle(none)color(none)) nobox )leg(on order(4 1 5 2 6 3))xtitle("Month") ytitle("log(UNITS)") ylabel(-2(1)2))(scatter coef month if type=="s" & program=="hu15" & label=="a2", connect(l) lpattern(solid) lcolor(red) msize(small)  mcolor(red)  msym(O) xlabel(-3 "Jul."-2 "Aug." -1 "Sep." 0"Oct." 1 "Nov." 2 "Dec." 3"Jan." 4 "Feb." 5 "Mar." 6 "Apr.")) (scatter coef month if type=="s" & program=="hu15" & label=="a1", connect(l) msymbol(T) msize(small))

graph export "HU_15_label2.eps", as(eps) fontface("Times New Roman") replace


***Subsidized prices by labels A+++, A++, and A+, time path, window around reform [-3,3].
twoway(rarea max95 min95 month if type=="sp" & program=="hu15" & label=="a3", color(gray%15)) (rarea max95 min95 month if type=="sp" & program=="hu15" & label=="a2", color(green%50) lcolor(green%50) lpattern(dash)) (rarea max95 min95 month if type=="sp" & program=="hu15" & label=="a1", color(none) lcolor(black) lpattern(dot)) (scatter coef month if type=="sp" & program=="hu15" & label=="a3", connect(l) tline(0, lstyle(refline) lcolor(black) lwidth(0.1)) tline(3, lstyle(refline) lwidth(0.1) lcolor(black) lpattern(dash)) yline(0, lstyle(refline) lcolor(gs8) lwidth(0.1)) lpattern(solid) lcolor(black) mcolor(black) msym(D) msize(small) graphregion(color(white)) legend(ring(0) position(10) row(3)  lab(4 "A+++")lab(1 "95% CI")lab(5 "A++") lab(2 "95% CI") lab(6 "A+")lab(3 "95% CI") region(lstyle(none)color(none)) nobox )leg(on order(4 1 5 2 6 3))xtitle("Month") ytitle("log(PRICE)") ylabel(-0.05(.05)0.1))(scatter coef month if type=="sp" & program=="hu15" & label=="a2", connect(l) lpattern(solid) lcolor(red) msize(small)  mcolor(red)  msym(O) xlabel(-3 "Jul."-2 "Aug." -1 "Sep." 0"Oct." 1 "Nov." 2 "Dec." 3"Jan." 4 "Feb." 5 "Mar." 6 "Apr.")) (scatter coef month if type=="sp" & program=="hu15" & label=="a1", connect(l) msymbol(T) msize(small))

graph export "HU_15_label_p2.eps", as(eps) fontface("Times New Roman") replace





